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Abstract. We consider the matrix 3p = Zp + Zp, where the entries of Zp are the 
values of the zeta function of the finite poset P. We give a combinatorial interpreta- 
tion of the determinant of 3p and establish a recursive formula for this determinant 
in the case in which P is a boolean algebra. 



II. Introduction 



The theory of partially ordered sets (posets) plays an important role in enumer- 
ative combinatorics and the Mobius inversion formula for posets generalizes several 
fundamental theorems including the number-theoretic Mobius inversion theorem. 
For a detailed review of posets and Mobius inversion we refer the reader to [SI], 
chapter 3, and [Sa]. Below we provide a short exposition of the basic facts on the 
subject following [SI]. 



A partially ordered set {poset) P is a set which, by abuse of notation, we also call 
P together with a binary relation, called a partial order and denoted <, satisfying: 

(1) x < x for all x G P (reflexivity). 

(2) If x < y and y < x, then x = y (antisymmetry). 

(3) If x < y and y < z, then x < z (transitivity). 
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Two elements x and y are comparable if x < y or y < x. Otherwise they are 
incomparable. We write x < y to mean x < y and x 7^ y. 

Examples. 

(1) Let nGN. The set [n] = {1,2,... , n} with the usual order forms a poset 
in which any two elements are comparable. Such a poset is called a chain 
or a totally ordered set. 

(2) Let n G N. Consider the poset P n of subsets of [n] under the inclusion 
relation. This poset is called a boolean algebra of rank n. In [SI] it is 
denoted by 2H 

(3) Let n G N. The set _D n of all positive divisors of n forms a poset under the 
order defined by i < j in D n if j is divisible by 2. This poset is called the 
divisor poset. 

A closed interval [x, y] is defined whenever x < y by [x,y] = {z <E P : x < z < y}. 
The empty set is not regarded as an interval. Denote by Int(P) the set of intervals 
of P. If / : Int(P) -> C, we write for 

The incidence algebra I(P) of P is the C-algebra of all functions / : Int(P) — > C 
under the convolution given by 



The zeta function £ G / (P) of a poset P is defined by ((x, y) = 1 for all x < j/ in 
P. If P is a locally finite poset (i.e. every interval in P is finite), the zeta function 
C is invertible in the algebra I(P). Its inverse is called the Mobius function of P 
and is denoted by [i. Note that one can define fx inductively by 

fj,(x,x) = 1, for all x G P, 



Proposition 1. (Mobius Inversion Formula) Let P be a finite poset and 




x<z<y 



It is known that I(P) is an associative C-algebra with identity 





fj,(x, z), for all x < y in P. 



x<z<y 



f,g:P^C. Then 



9(x) = ^2 f(y) f° r aUxe f(x) = ^ 9(y)Ky, x) for all x G P. 

y<x y<x 
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Proof. See [S], pg. 116. □ 

There is also a dual form of the Mobius Inversion Formula. Moreover, the formula 
works for more general posets. All that is needed is that every principal order ideal 
be finite. 

Remark. If P is the divisor poset of Example (3) , the proposition above becomes 
the number-theoretic Mobius Inversion Theorem. If P is the boolean algebra of 
Example (2), the proposition gives the Principle of Inclusion-Exclusion. Finally, if 
P is the chain of Example (1), the proposition becomes the Fundamental Theorem 
of the Difference Calculus. For the aforementioned theorems see [Sa] . 

For the remainder of the article, P will be a poset with n elements and the 
partial order denoted by <. We choose a labelling x\, X2, . . . ,x n of the elements of 
P such that Xi < Xj =>- i < j. 




Definition. The zeta matrix Zp of a poset P is defined as the n x n matrix with 
entries 



Observe that, with the chosen labelling, the zeta matrix is unipotent upper 
triangular. Its non-zero entries are the values of the zeta function. 

We define the matrix 3p by 3p = Zp + Z P . In Section §2 we give a combinatorial 
interpretation of the determinant of 3p- Section §3 refines the interpretation in the 
case in which P is a chain. The value of the determinant in this case is n + 1. The 
main theorem of the paper evaluates the determinant of 3n := 3p„ when P n is the 
boolean algebra of rank n. More specifically, in Section §4, we prove the following 
recursive formula on n. 

Main Theorem. If n > 3 is odd, then det(3 n ) = 0. If n is even, then det(3 n ) = 
2 Q ™ where «2 = 2, and a n = 4ct n _ 2 — 2 for n > 4 and even. 

Consider also the matrix 97lp defined by Tip = Mp + M P , where Mp = Zp 1 . 
The non-zero entries of M P are the values of the Mobius function. We refer to M P 
as the Mobius matrix of the poset P. We have the following theorem. 

Theorem 1. det(97t P ) = det(3p). 

Proof. The theorem is a direct consequence of the following lemma. □ 
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Lemma 1. Let U be an n x n matrix such that det(U) = 1 and let V = U 



Then, det(t/ + IP) = det(V + V f ). 
Proof. We have 

v* + v = (u- 1 ) 1 + u- 1 = (u-yuu- 1 + (u*)- 1 ^- 1 = (u-yp + iPp- 1 . 

Thus det(V + V*) = det((t/" 1 ) t ) det(U + U*) det(t/" 1 ) and since det(*7) = 1, we 
have det(y + V*) = det(t/ + t/ t ). □ 



In this section, we will give two related combinatorial interpretations of det(3p)- 
The first involves collections of disjoint cycles in a directed graph Dp associated to 
P. The second involves cycle decompositions of elements of a certain subgroup S% 
of the symmetric group S n associated to P. 

Consider a poset P as in the previous section with \P\ = n. The matrix Yp = 
Zp — I n , in which the diagonal entries of Zp are replaced by 0, can be interpreted 
as the adjacency matrix of a directed graph (digraph) Gp associated to the strict 
order relation x < y in P. The vertices of Gp are the elements of P, and there is a 
directed edge from x to y if and only if x < y. Similarly, the matrix %)p = Yp + Y P 
is the adjacency matrix of the directed graph Dp in which there are edges in both 
directions between each pair of distinct comparable elements x,y G P. Then we 
have 



and det(3p) = det(2)p + 2I n ) = %(— 2) where x(t) is the characteristic polynomial 
x(t) = det(2)p — tl n ) of the matrix 2)p. As in [M], we will also call this the 
characteristic polynomial of the graph Dp. 

If P = {xi, . . . , x n } and 2)p — tl n = (m^), then by the definition of the deter- 
minant 



§2- 



Combinatorial interpretations of det(3p) 



3p = Z P + Z P = Yp + Y P + 2J n = fQ P + 21, 



n 



x(t) 



E( 



Then we have 



n 







if Xj not comparable to x a ^ for some j, 
otherwise, 



n m Mj) 



3=1 
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where k a is the number of fixed points of the permutation a. 

Let S£ = {a E S n : xj comparable in P to x a (j), for all j = 1, . . . , n}. Then 
combining the last two displayed formulas 

(i) x (t)= (-iy i9n ^\-t) k °. 

As shown in the work of Harary and Mowshowitz (see [H] and [M]), there is an 
interesting relation between the coefficients of %(t) and cycles in the graph Dp. 
Namely, if we write 



X 



i=0 



then ao = 1 (the coefficient of (— t) n equals 1 since the identity permutation in 
S n is the only permutation with exactly n fixed points) and Theorem 1 of [M] 
demonstrates that for i > 1, 




(2) * = i s K-ir j+1 1 f DP (iu...,ir) 



where the summation is taken over all partitions i = %\ + • • • + i r with ij > for 
all j, and fD P (h, ■ ■ ■ , ir) is the number of collections of disjoint directed cycles in 
L>p of lengths ii, . . . , i r . Thus 



H = ^(-1) 1+7 7z)p(h,--- ,ir)- 



Note that in fact ij > 2 for all j since Dp contains no loops (no edges from a 
vertex to itself). Therefore a± = 0, which corresponds to the fact that there are no 
permutations in S n with exactly n — 1 fixed points. 

Example. To illustrate the relation between (1) and (2), let P = P 2 be the boolean 
algebra of rank 2. The ordering of the elements of P is given by the correspondence 

0< — >00, {1} < — ►01, {2} < — ► 10, {1,2}< — > 11. 

Thus we have the labelling 

p = { Xl = 00, x 2 = 01, x 3 = 10, x 4 = 11}. 

Then 

(2 1 1 1\ 
12 1 
10 2 1 
Vl 1 1 2/ 
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and 



X (t) = det(2) P - th) = det 



f-t 1 1 l\ 

1 — * 1 

1 — t 1 

V i i i -*/ 



= ^a,(-t) 4 - = (-t) 4 -5(-t) 2 + 4(-t), 



i=0 



where the a^'s are defined as above, so ao = 1, ai = 0, ai = —5, a3 = 4, a4 = 0. 
The digraphs Gp and Dp are given below. 



G P : 




Dp: 




We obtain one term of (— t) 2 for each permutation in with exactly 2 fixed 
points. Since each part of the partition must be at least 2, the only allowable 
partition of 2 is 2 = 2 and fp> p (2) counts the following collections of cycles in Dp 
which correspond to the transpositions in S[: 



01 



11 




00 



•10 



01 



11 




o # o 



no 



11 



01« 




10 



00 



01« 



11 




# 



10 



01« 



11 

A 

V 

00 



•10 



This corresponds to the coefficient a,2 = 



-5 
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We obtain one term of (— t) 1 for each permutation in with exactly 1 fixed 
point. The only allowable partition of 3 is 3 = 3 and fo P (3) counts the following 
collections of cycles in Dp which correspond to the 3-cycles in S[: 



11 11 




This corresponds to the coefficient as = 4. 

Finally, we obtain one term of (— 1)° = 1 for each permutation in with no 
fixed points. Such permutations are either 4-cycles or products of 2 transpositions. 
The only allowable partitions are 4 = 4 and 4 = 2 + 2. Then /d j ,(4) counts the 
following collections of cycles in Dp (with a — sign) which correspond to the 4-cycles 
in Si': 



11 11 




00 00 



On the other hand, /p, p (2, 2) counts the following collections of cycles in Dp 
(with a + sign) which correspond to the permutations in S± which are products of 
2 transpositions: 



11 11 




8 



DETERMINANTS OF ZETA MATRICES 



Thus a 4 = -2 + 2 = 0. 



Returning to our general situation, substituting from (1) we have 

n 

X (-2)=5>2— 

= 2n + E( E (-l) i+r fD F (iu...,ir))2 n - i . 

4=1 \E*i=< / 

For each collection of r disjoint directed cycles there are n — i vertices not contained 
in the cycles. Let s = n — i + rbe the total number of connected components of 
the union of the cycles and the disconnected vertices. We reindex the last sum as 
a sum over s: 



( 



x(-2) = $>iy 



s=l 



\ 



E fopiii, ■ ■ ■ , v)2 4 



,(;,■) 



y s = r-\~n — 



E(- 1 ) n " 



s=l 



(The sign comes from the fact that i + r = s — n + 2i, so (— l) l+r = (— l) n s .) The 
coefficient 

Cs= fDp(h,... Jr)2 S ~ r 

r,(ij) 
s = r-\-n — ^2 ij 

can also be interpreted as the integer 



(where there are s — r ones) counting the number of collections of disjoint directed 
cycles of lengths i\, . . . , i r , 1, . . . ,1 containing all the vertices in a directed graph 
Dp obtained from Dp by adding two loops at each vertex (corresponding to the 
diagonal entries 2 in the matrix 3p = Zp + Z P ). Hence c s is the number of 
collections of disjoint directed cycles in the graph Dp having precisely s connected 
components and containing all the vertices. Note that s — r is the number of loops 
contained in each such collection. 

To summarize, we have proven the following theorem. 



DETERMINANTS OF ZETA MATRICES 



9 



Theorem 2. With the above notation, for all finite partially ordered sets P, 

n 

det(3p) = ^(-l) n " s c s . 



s=l 



In terms of permutations in S% , ifl<s<n — 1, we have 



(3) 



and c n = T 



8-1 

= E 



a g g p . °" i s a product of s — £ disjoint 
n ' cycles of length > 1 



§3. The chain case 



Consider the special case in which P is the chain of length n. Then P = [n] = 
{1,2,... , n} with the usual ordering. All elements of P are comparable, and we 
have 



Z\ n ] = 



(1 1 ••• 1\ 
1 ••• 1 



1/ 



(1 -1 ••• \ 
1 -1 ... 



and M [n] = = 



\ 



1 -1 
1 / 



where all lower-triangular entries in each matrix are zero. 



Lemma 2. With the above notation 



det(3[„]) = det(9Jl [n] ) = Ti+l. 



Proof. The first equality holds by Theorem 1. We prove the second equality by 
induction on n. We have 9DT[i] = (2), hence the lemma holds when n = 1. Suppose 
the lemma holds for all k < n, and consider 



2%1 



/ 2 -1 
-1 2 -1 

-1 '•• 



\ 



2 -1 
-1 2 / 
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where all remaining entries are zero. Expanding along the first row, we have 
det(97l[ n ]) = 2det(9Jl[ n _i]) - (-1) det(B n ) where 



B n = 



/- 1 - 1 


V 



... o 



Wl 



[n-2] 



J 



Now expanding along the first column in B n , we have det(B n ) = — det(9Jt[ n _ 2 ]). 
Thus we have 

det(SD r t[ n ]) = 2 det(9%_i]) - det(9Jl [n _ 2] ) = 2n - (n - 1) = n + 1 
by applying the inductive hypothesis. □ 



All elements of P are comparable in the chain case, and so the digraph Dp 
is complete with double edges between all vertices and double loops at all ver- 
tices. Considering permutations written as a product of disjoint cycles, denote 
by r(n; 71, 72, • • • , 7n) the number of permutations of [n] with exactly 7$ cycles of 
length i for each 1 < i < n. By [SI], Proposition 1.3.2 we have 



T(n; 71,72,... ,7 n ) 



1^71! 2^72! ... n^ ln V 



Then, in the chain case, the formula for C{ in (3) becomes 
Ci= r (^;7i,72,--- ,7n)2 71 . 

(71 , ■ •• ,7n) 
E 7,"=i 



The Cj's may be interpreted as modified Stirling numbers of the second kind. 

count the number of 



Recall that the Stirling numbers of the second kind 



m 



permutations of n elements with m disjoint cycles and are given by 



n 
m 



y 

^ l mi mi! 

E 

E im i 



■m,; =m 



mi! 2 m2 m 2 ! . . . k mk m k l 



Y T(n;mi, . . . ,m k ). 



E m i= m 

E im i= n 



Corollary 1. The alternating sum of the modified Stirling numbers satisfies the 
following. 



£(_!)»-* r(n; 7 i,...,7»)2 7l =n+l 



i=i 



E 
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§4. The Boolean algebra case 



Let [n] = {1, 2, . . . , n} and consider the poset P n = 2^ of subsets of [n] under 
the inclusion relation. Let 3™ be the matrix defined in §1. The main result of this 
section if the following theorem. 

Theorem 3. If n > 3 is odd, then det(3n) =0. If n is even, then det(3n) = 2 a ™ 
where = 2, and a n = 4a n _ 2 — 2 for n > 4 and even. 



Proof. The proof will rely on several lemmas. First, we identify a particularly 
useful labelling of P n = 2^ for our purposes, and we consider only this labelling 
in the following. Each subset A e P n will be encoded as a binary vector v(A) of 
length n : 

v(A) = (v n ,v n -i,... ,vi) 

where 

f 1 ifieA, 
I if not. 

Our labelling of P n induces the usual numerical ordering when we interpret v(A) 
as the binary expansion of an integer m, with < m < 2 n — 1. 

Using this labelling yields an interesting recursive structure in the matrices 3n- 

Lemma 3. The matrices Z n and 3n have the following properties. 

(1) The entries of Z n above the diagonal are the first 2 n rows in the Pascal 
triangle modulo 2. 

(2) For n > 2, Z n and 3n have block decompositions: 

Zn—l Z n —\ \ , _ ( 3n-l Z n -\ 

and i). 



Z n -i J " \Zn-l 3n-l 

(This statement also holds with n = 1 if we take 3o = 2, Z = 1.) 
(3) The Z n matrix sequence can be generated by a recursive procedure as follows. 

Given Z n _\, to form Z n we replace each entry 1 by a 2 x 2 block 



1 

and each entry by a 2 x 2 zero matrix. From part (1), we get a similar 
recursive procedure for the 3n sequence. 



Proof. These properties follow directly from the definition of the 3n matrices and 
the properties of the preferred ordering on P n . □ 

To evaluate the determinant det(3n), we follow the general advice of [K] and 
introduce parameters in the matrix entries. Specifically, we consider the matrix: 
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Using the Maple computer algebra system, the determinants of the first few of these 
matrices are found to be: 



n 



det(3 n (x,y)) 



(4) 



1 

2 
3 
4 
5 



x + xy + y 
(x + y) 2 (x 2 — xy + y 2 ) 
(x - y) 2 (x 2 + + y 2 ) 3 
(x + y) 6 (a; 2 - xy + y 2 f 
(x-y) 10 (x 2 + xy + y 2 ) 11 



An interesting recurrence explains the patterns evident in these examples. 

Lemma 4. The determinants of the "5 n ( x ,y) matrices are related by the following 
recurrence: 

det(3n+2(>,2/)) = (det(3n(£,?/)) 2 det(3n+i(>, -y)). 

(Note that the right side involves both 3n+i and 3n and that y is negated in the 
second factor.) 



/xZ n + yZ l n 


xZ n 


xZ n 


xZ n 




xZ n + yZ l n 





xZ n 


yzi 





xZ n + yZ^ 


xZ n 


V yzi 


yzi 


yzi 


xZ n + yZ l n 



Proof. We use part (2) of Lemma 3 twice to form the following block decompo- 
sition of 3 n+2 (x,y): 



3 n +2(x,y) = 



(where each entry is a block of size 2 n x 2 n ). To evaluate the determinant, we 
perform block-wise row and column operations. To simplify the notation, we write 
Z = Z n First subtract row 4 from each of the first three rows to obtain: 

( xZ xZ - yZ l xZ - yZ l -yZ l \ 
xZ -yZ l -yZ l 

-yZ* xZ -yZ l 



\yz l 



yz l 



yz l 



X 



Z + yZ 1 / 



Subtract column 1 from each of the columns 2,3,4 to obtain 

-xZ-yZ*\ 

-yz l 
-yz l 

xZ 



( xZ -yZ l -yZ l 

xZ -yZ l 

-yZ* xZ 

\yZ l 



Next, subtract column 2 from column 4 



/ xZ -yZ l -yZ f 

xZ -yZ l 

-yZ* xZ 

\yZ l 



-xZ \ 

—xZ — yZ l 
' 
xZ 



) 
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Add column 1 to column 4, then add row 4 in the resulting matrix to row 2: 

\ 



3' = 



( xZ -yZ 1 -yZ 1 

yZ 1 xZ -yZ 1 

-yZ* xZ 

\yZ l xZ + yZ 1 ) 



Expanding along the last column we obtain 

xZ -yZ 1 

(5) det(3n+ 2 (£, y)) = det(3') = det(xZ + yZ f ) det | yZ l xZ 

-yZ 1 



-yz 1 ' 

xZ 



The first factor on the right of (5) is det(^> n (x,y)). Continuing with the 3x3 
matrix, subtract column 3 from column 2: 

xZ -yZ* ' 

yZ l xZ + yZ l -yZ l 

— xZ — yZ l xZ 

then add row three to row 2: 

xZ -yZ* 

yZ f xZ- yZ f 

-xZ - yZ f xZ 



Expanding along column 2, we have 

det(3n+2(>, y)) = det(3n(>, y)) 2 det 

= det(3n(^,y)) 2 det 



xZ —yZ 1 

yZ 1 xZ - yZ 1 

xZ - yZ 1 -yZ 1 
xZ xZ — yZ 1 



= det(3n(>, y)) 2 det(3n+i(£, -y)), 

as claimed. For the last equality, we perform row and column interchanges to put 
the final matrix shown into the form: 

xZ — yZ 1 xZ 
-yZ 1 xZ - yZ 1 

required for 3n+i(x, —y). □ 

From the initial cases computed with Maple in (4) and the recurrence from 
Lemma 4, we see that there are nonnegative integers a n , (3 n such that 

(6) det(3„(x, y)) = (x + (-l) n y) a -(x 2 - (-l) n xy + y 2 f- . 

Moreover, the recurrence from Lemma 4 implies that 

^ f an+2 = 2a n + a n+1 , 

I Pn+2 = 2/3 n + Pn+l- 

We also have the following fact that is evident from (4): 
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Lemma 5. For all n > 1, a n = f3 n + (— l) n+1 . 

Proof. This follows by induction. The base cases come from the Maple com- 
putations in (4) above: ol\ = 0,/3i = 1, and a 2 = 2, /?2 = 1- For the induction 
step, assume that the claim of the lemma has been proved for all I < k + 1. Then 
subtracting the two recurrences from (7) shows that 

a k+2 -(3 k+2 = 2(a k -p k )+a k+1 -p k+1 = 2(-l) k+1 + (-l) k + 2 = (-l) k+1 = (-l) fc + 3 . 

We are now ready to conclude the proof of our main theorem. To determine 
the determinant of the original 3n = 3n(l> 1)) we simply substitute x = y = 1 
in (6). The factors of x — y show immediately that det(3n) = if n is odd. 
Moreover, when n is even we have det(3 n ) = 2 a ™. We solve the first recurrence 
in (7) for a n by the standard method for second order linear recurrences with 
constant coefficients. The characteristic equation is r 2 — r — 2 = 0, whose roots 
are r = 2, —1. Hence a n = ci(2) n + ci{— l) n for some constants ci, C2- The initial 
conditions a.\ = 0, a 2 = 2 show that c\ = 1/3, c 2 = 2/3. Hence: 

«» = ^(2) B +|(-1) B . 
Hence if n, and therefore also n + 2, are even, we have 

a n+2 = ^(2r+ 2 + | = 4 Q(2)» + |) - 2 = Aa n - 2. □ 



Corollary 2. 7/n even, then 

2™+2 

det(3n) = 2^~. 
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